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Abstract 

Ultrathin Bi 2 Se 3 -NbN bilayers comprise a simple proximity system of a topological insulator and 
an s-wave superconductor for studying gating effects on topological superconductors. Here we 
report on 3 nm thick NbN layers of weakly connected superconducting islands, overlayed with 
10 nm thick Bi 2 Se 3 film which facilitates enhanced proximity coupling between them. Resistance 
versus temperature of the most resistive bilayers shows insulating behavior but with signs of 
superconductivity. We measured the magnetoresistance (MR) of these bilayers versus temperature 
with and without a magnetic field H normal to the wafer (MR=[R(H)-R(0)]/{[R(H)+R(0)]/2}), 
and under three electric gate-fields of 0 and ±2 MV/cm. The MR results showed a complex set of 
gate sensitive peaks which extended up to about 30 K. The results are discussed in terms of vortex 
physics, and the origin of the different MR peaks is identified and attributed to flux-flow MR in 
the isolated NbN islands and the different proximity regions in the Bi 2 Se 3 cap-layer. The dominant 
MR peak was found to be consistent with enhanced proximity induced superconductivity in the 
topological edge currents regions. The high temperature MR data suggest a possible pseudogap 
phase or a highly extended fluctuation regime. 

Keywords: topological superconductivity, proximity effects, gating effects, topological insulator - 
superconductor bilayers 


1. INTRODUCTION 


Surface edge states of topological superconductors 
(TOS) are expected to support zero energy modes or 
Majorana fermions which are robust against disorder 
and decoherence mm- These zero energy excitations 
could therefore be useful in potential application in 
spintronics and quantum computing [Sill!. Bulk TOS 
such as copper doped Bi 2 Sez should have been the 
simplest materials to study TOS properties, but compli¬ 
cations due to their inherent inhomogeneity |5| and the 
presence of possible superconducting impurity phases 
such as CuS 62 |5|, make them less attractive for such 
investigations. An alternative way for realizing TOS is 
by inducing superconductivity in a topological insulator 
or in semiconductor-nanowires with strong spin-orbit 
interaction via the proximity effect (PE) [TMO]. Uncon¬ 
ventional superconductivity in these systems, such as 
revealed by the presence of zero bias conductance peaks 
(ZBCP), indicates zero energy bound states that might 
be due to Majorana zero energy modes, but could also 


originate in zero energy Andreev bound states. It is hard 
to distinguish between these two different phenomena, 
and efforts are ongoing in order to achieve this goal 

mma. 

Spatial sharpness of the boundary region between the 
superconductor and the topological or semi-conducting 
material is also essential in order to distinguish between 
the near-zero-energy end states originating in Andreev 
bound states, and the Majorana zero energy modes 
in the topological case m- Since this boundary is 
generally created by gating in the experiments, its un¬ 
avoidable gradual spatial change adds more uncertainty 
to the interpretation of the observed ZBCPs as due to 
the Majorana modes imniiii]. The role of gating in 
these nanowire-superconductor experiments was further 
investigated and a variety of additional phenomena such 
as ZBCP oscillations versus gate voltage and magnetic 
field were observed and interpreted in the context of 
Majorana modes as well as alternatives such as Kondo 
and disorder effects [1^]. Hence, gating is sufficiently 
important in these studies and we decided to use a 
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FIG. 1: (Color online) Atomic force microscope (AFM) im¬ 
age of nominally 3 nm thick NbN layer, together with a height 
profile along the line shown in the image. The two horizon¬ 
tal lines near the 9-2=7 nm tall particulate are artifacts due 
to the AFM inability to follow smooth traces near a big dis¬ 
turbance (in tapping-mode here while scanning along the x- 
direction). 





FIG. 2: (Color online) AFM image of a 10 nm 812863 on a 
3 nm NbN bilayer (a), together with two height profiles (b) 
and (c) along the long and short lines shown in the image, 
respectively. 


simple proximity system of a bilayer comprising of a 
topological insulator and a 2 D superconductor for study¬ 
ing gating effects on topological superconductors [TB] . In 
continuation to our previous studies of 8128— NbN 
junctions [m US] we used ultra-thin bilayers of this 
system, even thinner than used before m, where 
the superconducting NbN islands with weak-links in 
between them are overlayed with a thicker Bi2Se3 layer 
which enhances the coupling between the NbN islands 
via the (inverse) proximity effect. Here we report on 
gating effects on the magnetoresistance (MR) of this 
hybrid system, which shows a highly gate sensitive, 
non-monotonous behavior versus temperature. 


2. EXPERIMENTAL METHODS 


2.1. Preparation and surface morphology of the 
films and bilayers 


The NbN and 812863 thin films were prepared 
as described in detail previously m- Briefly, laser 
ablation deposition was used where the NbN films were 
deposited under 30 mTorr of N2 gas flow and at 600 
°C heater block temperature, while the 812863 layers 
were deposited under vacuum and at 300 °C. All films 
and bilayers were deposited on fused silica wafers and 
yielded grainy films with rms roughness of ~ 10 % of the 
films thickness. Fig. 1 shows an atomic force microscope 
(AFM) image of the surface morphology of a 3 nm thick 
NbN film, together with a typical height profile along 
the line shown in the image. While the lateral size of 
the NbN grains is ~ 50-100 nm, one can see larger bright 
areas of at least 1000 nm in size which represent thicker 
superconducting NbN islands. These are separated by 
broad darker areas which constitute thinner weak-links 
between these islands. This islands structure that 
comprises a network of strong superconducting regions 
connected by weak-links is essential for the present 
study, otherwise a superconducting short will mask all 
our transport data. SEM images of similar NbN films 
on glass showing their grainy nature can also be seen in 
m- In addition. Fig. 1 shows a few 1-3 nm deep holes 
of about 60 nm diameter while the overall rms roughness 
of this film is ^ 0.35 nm. 


All the 812863 — NbN bilayers in the present study 
were obtained using an ’’in-situ” process where both 
the 3 nm thick NbN film and the 10 nm thick 812863 
cap layer were prepared in the same deposition run 
without breaking the vacuum. This kept the interface 
between the NbN and 8i28e3 layers protected against 
contamination and oxidation which occur if the NbN 
surface is exposed to air [THUD]. Fig. 2 shows an AFM 
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image of the surface morphology of such a bilayer (a), 
together with two height profiles (b) and (c) along the 
long and short lines shown in (a), respectively. The 
bilayer shows taller grains on a smoother background, 
with grains height of 4-8 nm and some holes, the deepest 
of which is of about 5 nm deep. Shallower holes are 
also visible which reveal a step-like depth profile as seen 
in (c), where the step heights are of about 1 or 2 nm 
in agreement with one or two quintuples of the Bi 2 Se 3 
structure. The overall rms roughness of this bilayer is 
^1.1 nm. 

2.2. The gates and contacts geometry 


Transport measurements were done using an array 
of 40 gold coated spring loaded spherical tips for the 
4-probe measurements on 10 different locations on the 
wafer (Ci with i=l to 10, 4 contacts for each location). 
Fig. 3 depicts a schematic drawing of the bilayer, gates 
positions and 20 representative contacts. Contacts Cl 
and CIO are separated from the bilayer by two scratches, 
while the 4 contacts of each are shorted and connected 
to the top and bottom gates, respectively. Since the 
fused silica substrate is 1.5 mm thick and its dielectric 
constant is low, the bottom gate is useless here and is 
shorted to the bilayer and ground, as can be seen in 
the schematic cross-section of the bilayer and gates in 
Fig. 4. The insulation layer of the top gate is made of 
a much thinner PMMA resist of 1.3 p,m thickness. This 
allows for an electric field of ±2 MV/cm to be obtained 
when a voltage of =Fl00 V is applied to the top gate 
(opposite polarity). This will be marked in the following 
by =f 100 Vg. Fig. 4 shows two weakly connected NbN 
islands (red dots) together with their depletion layer 
(blue dots on white background) in the Bi 2 Se 3 layer 
near the interface with the NbN islands. Proximity 
induced superconductivity also occurs in this layer and 
can extend up to the top surface of the Bi 2 Se 3 film m- 



FIG. 3: (Color online) (a) A schematic drawing of the Bi 2 Se 3 
on NbN bilayer, together with the top and bottom gates and 
some of the Ci contacts where i=l to 10 (each Ci has 4- 
point contacts). The Cl and CIO contacts are separated from 
the bilayer and dedicated for contacting the top and bottom 
gates, respectively. 
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FIG. 4: (Color online) A schematic cross-section of the bi¬ 
layer and gates. Two superconducting NbN islands marked 
by red dots are weakly connected by the non-superconducting 
green shaded area. Helical edge currents for the two spin ori¬ 
entations of the Bi 2 Se 3 layer are depicted by the horizontal 
arrows on the top and bottom surfaces. The depletion layer 
in the Bi 2 Se 3 is marked by the blue-dots. The back gate 
is grounded and shorted to the bilayer in the present study, 
and a negative top gate voltage -Vg corresponds to a positive 
electric field E. 
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3. RESULTS AND DISCUSSION 
3.1. Comparing bilayers to NbN films 

We start this study by repeating the main experiment 
of our previous work [16] of comparing the properties of 
a Bi 2 Se 3 -NbN bilayer to those of a reference NbN film. 
This time however, half as thin bilayer and reference 
film are used, in order to enhance the weak-link behavior 
of both and minimize superconducting shorts. No gates 
were prepared on the wafer used for this experiment in 
order to minimized its exposure to ambient air which 
is detrimental to the NbN film [20]. A bilayer of 10 
nm 812863 on 3 nm NbN was deposited on half the 
fused silica (FS) wafer, while the 3 nm thick NbN 
reference film was deposited on the other half. The 
two halves were separated by a scratch as shown in 
Fig. 2 (a) of Ref. [T^. Scratch separation rather 
than photolithography was used in order to avoid 
deterioration of the films. The transport results on 
this wafer are presented in Figs. 5 and 6. As observed 
previously, the R versus T results show that the onset 
Tc of the reference NbN film is higher by about 1 K 
than that of the bilayer. This indicates a standard 
proximity effect (PE) where the normal electrons of the 
B 12863 suppress superconductivity in the NbN layer 
of the bilayer [16]. The inverse PE however, where 
enhanced superconductivity in the bilayer occurs due to 
pair penetration into the 812863 , is absent now, at least 
down to the minimum temperature used in this study 
(1.8 K). This, together with the fact that both bilayer 
and reference film do not reach zero resistance down to 
1.8 K, is a good indication that the weak-links in both 
are actually weaker, as originally planned. The inset 
to Eig. 5 shows the magnetoresistance of the bilayer 
and reference film versus temperature. It shows that 
the Tc onset (determined by the rise of the MR) of the 
reference film is again, larger than that of the bilayer. 
Fig. 6 exhibits a zoom-in on the low MR data and its 
extension to higher temperatures. It shows that the MR 



FIG. 5: (Color online) Resistance versus temperature of a 
pristine, half bilayer (BL) - half reference layer (RL) on a 
fused silica wafer under 0 Vg and 0 T. The BL is a 10 nm 
Bi 2 Se 3 on 3 nm NbN bilayer and the RL is a 3 nm NbN film. 
Also shown is RL data normalized to the BL data at 10 K 
(norm C7 RL). Inset: The magnetoresistance at 0 and 1 T 
magnetic field normal to the wafer versus temperature. 

onset of the reference NbN layer (C7 RL) is at about 
5 K, in agreement with the R versus T data of Fig. 5. 
It also shows that above 5 K, its MR is zero to within 
the noise of the measurement. This is different from the 
bilayer data of Fig. 6 (C2 and C4), where a non zero, 
positive MR signal of a few percents is observed up to 
50 K. For comparison we show an MR signal of a 10 nm 
thick 812863 him on FS, which is very similar to the 
bilayer data. We thus conclude that above 5 K in this 
bilayer, the MR signal originates in the 812863 cap layer 
with no noticeable effect of the underlying NbN layer. 

The inset to Fig. 6 shows the resistance versus 
magnetic held H normal to the wafer at 1.88 K, for 
the same bilayer and reference him contacts. All the 
resistances are increasing with held, but while the bilayer 
resistance tends to saturate at high helds, that of the 
reference him increases quite linearly with H. The later 
indicates hux-how resistance in a superconductor as the 
mechanism for the linearly increasing resistance. If we 
go back now to the inset of Fig. 5, we see that this 
effect is stronger in the reference NbN him than in the 
bilayer. The peak MR at ~1.9 K in all contacts is due 
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FIG. 6: (Color online) Main panel: Zoom in on the high 
temperature tail of the MR data of the bilayer and reference 
film of the inset to Fig. 5. Also added for comparison is the 
MR versus T of a bare 10 nm thick film of Bi 2 Se 3 on fused 
silica. Inset: Resistance versus magnetic held at 1.88 K of the 
same BL and RL contacts. 

to maximum flux flow which results from compensation 
between increased vortex generation and pinning effect 
with decreasing temperature. Another interesting result 
in the inset to Fig. 5, is that in the C2 BL contact 
there is a clear knee in the MR data at T~2.1 K which 
indicates two different MR behaviors below and above 
it. We attribute the MR above this temperature to 
the stronger superconducting NbN islands as seen in 
Fig. 1, while below it, where the weak-links (the darker 
areas in Fig. 1) also become superconducting due to 
proximity coupling via the Bi 2 Se^, the MR increases 
further on lowering the temperature. In the following we 
shall see more supporting evidence for this interpretation. 

3.2. R and MR of reference Bi 2 Se 3 films 

For understanding the MR of the bilayers one needs 
in addition to the MR of the 3 nm NbN reference film, 
also the MR behavior of a second reference film of 10 nm 
Bi 2 Se 3 . Such films were prepared and Fig. 7 presents 
resistance versus temperature results of a 10 nm Bi 2 Se 3 
film under 0 and 1 T magnetic field and -100 V gate 



FIG. 7: (Color online) Resistance versus temperature of a 
10 nm thick Bi 2 Se 3 reference film on FS with and without 
a magnetic held of IT. The inset shows the phase coherence 
length calculated for another 10 nm Bi 2 Se 3 him assuming 
weak anti-localization (WAL) and using the HLN model [24] . 

voltage. The resistivity of this film is ~6 mflcm or Okfl 
per square, which corresponds to an electron density 
of about 10^^ cm’ -3 EH. Thus the electron doping of 
this film is quite low, which results from the Se rich 
target used in its deposition process [16]. The R versus 
T data of Fig. 7 was used in order to calculate the 
magnetoresistance of this film and the result has already 
been presented in the main panel of Fig. 6. We found 
that MR versus T curves measured on the Bi 2 Se 3 film 
were insensitive to the gate voltage used here (tIOO Vg 
or ±2 MV/cm) to within the noise of the measurements. 
Steinberg et al. [22] did observe significant gating effects 
in a 20 nm thick Bi 2 Se 3 film, but their maximum 
electric field was about 40 times higher than what we 
used here. Thus, the presently used maximum gate 
voltage was insufficient to have any significant effect on 
the MR of the bare 10 nm Bi 2 Se 3 film. Fig. 8 shows the 
resistance versus magnetic field of this film at 1.85 K. 
The inset shows the data obtained on the pristine film, 
while the main panel shows the data of a one week old 
film which was kept under dry air. Clear aging effect was 
observed when the measurements were repeated after 
one week, as can be seen from the increased resistance 
of about 25%. Besides this, both sets of measurements 
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FIG. 8: (Color online) Resistance versus magnetic field at 
1.85 K for the 10 nm Bi 2 Se 3 film. The inset shows results 
measured on the pristine film while the main panel shows data 
obtained after the film was aged in dry air ambient for one 
week. Clearly, the resistance of the aged film went up due to 
deterioration of its surface layer. 


show a linear R versus H behavior at low fields up to 
about 0.1 T with saturation at higher fields. 



FIG. 9: (Color online) Resistance versus temperature of 

a reference bilayer of 10 nm An on 4 nm NbN under zero 
magnetic field and 0 Vg. 


can affect the resulting L^. Both however, decrease 
very rapidly versus T at low T and then decay more 
gradually. Again, this is not a proof that WAL occurs 
here, but if it does, these are the calculated ^^(T) curves 
from our measurements. 


The linear MR in Bi 2 Se^ is generally attributed to 
weak anti-localization (WAL) [^[^. If we assume that 
this is correct, we can use the HLN model [24] to cal¬ 
culate the phase coherence length as a function of 
temperature from our data. For a field of H = 1 T and 
in nm, the HLN model yields: 


1.2 X 10^ 


R{H) - i?(0) 
R{0)R{H) 


156, , ,156, 

^(2 + tt ) - ( 1 ) 

^ 4 , 


where tp is the digamma function. Since both R{H = IT) 
and R{0) were measured as a function of temperature, 
Eq. (1) allows us to extract L 4 which is thus also a 
function of temperature. The result is shown in the 
inset to Fig. 7 for another 10 nm thick Bi 2 Se^ film on 
FS. One can see that the resulting L 4 is very sensitive 
to the contact location on the wafer, probably because 
the coherent WAL scattering process is very sensitive to 
even slight inhomogeneities in the film. Also, since C2 
is closer to the edge of the wafer than C4, edge effects 


3.3. Reference bilayer of Au on NbN 


In this section we study the temperature dependence 
of R and MR of a 10 nm Au on 4 nm NbN bilayer on 
FS. This provides information on the proximity effect 
of NbN with a good metal, and allow for a comparison 
with the Bi 2 Se 3 -NhN bilayers which will be presented 
and discussed in detail in the next section. We chose to 
have a bit thicker NbN layer, of 4 nm instead of 3 nm, 
in this reference bilayer with 10 nm gold cap layer, since 
we worried that the ball-like gold grains will not fully 
protect the NbN layer from oxidation. The transport 
results of R and MR under different magnetic helds and 
gate voltages are shown in Figs. 9 and 10, respectively. 
The resistance increased gradually on cooling down from 
room temperature to a maximum at about 12-13 K, and 
its maximum value is about half that of the bilayer in 
Fig. 5 if it were prepared on a whole wafer. Fig. 9 
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shows a spread of the Tc values for the different contact 
locations on the wafer, with transition onset at about 6 
K. The MR data versus temperature of the 10 nm Au on 
4 nm NbN bilayer is shown in Fig. 10. The inset of this 
figure shows the MR at low temperatures under 2 and 
4 T magnetic fields, and under different gate voltages. 
The main panel is a zoom-in on the tail of the MR of 
the inset up to 90 K. The increase above about 1% of 
the MR in the inset marks Tc = 6 K of the main part 
of the bilayer which coincides with the resistive Tc onset 
value obtained in Fig. 9. Below this Tc value, the MR 
increases rapidly with decreasing temperature due to 
flux-flow in the bilayer. Between 4 and 5 K, one can 
observe a distinct knee in the MR which indicates the 
transition between superconductivity in the NbN islands 
above 5 K, and the proximity induced superconductivity 
in the gold layer in between islands below 4 K. The 
saturation MR at 200 % is a result of our definition 
of MR, which get this value when R(0)=0. These 
features are similar to those of the inset to Fig. 5 in 
a 10 nm Bi 2 Se 3 on 3 nm NbN bilayer under 0 Vg, 
where a lower Tc of 2.7 K was found, with a weaker 
signature of the NbN islands and inter-islands regimes, 
and without the saturation effect since the bilayer 
always remained resistive, even at 1.8 K. The inset to 
Fig. 10 also shows a small gating effect on Tc at 2 T, 
where under 0 Vg Tc is higher by ^0.3 K than at -100 Vg. 

The main panel of Fig. 10 shows that the MR above 
12-13 K and up to 90 K is within the noise of the 
measurements (±0.1 %). This is in contrast to the 
observation in the bilayer of Fig. 6 where an MR of 
a few percents is observed even much above Tc. If we 
determine Tc in Fig. 10 (somewhat arbitrarily) at T 
where MR> 0.2%, we find Tc values of 8 and 12 K at 
2 and 4 T helds, respectively. These Tc values mark 
the onset of flux-flow at the given fields, and are higher 
than the previously determined Tc from the inset (^6 
K). All these different Tc values reflect the fact that MR 



FIG. 10: (Color online) Magnetoresistance versus tempera¬ 
ture of the bilayer of Fig. 9 under different magnetic fields and 
gate voltages. The inset shows the low temperature behavior 
of the MR below Tc ~ 6 K. 


values higher than the noise level (reached at the onsets 
of flux-flow) occur at higher temperatures under higher 
fields, and depend sensitively on the onset criterion. 
In addition, the Tc = 12 K value coincides with the 
temperature value of maximum resistance as seen in 
Fig. 9. Therefore, we conclude that this is the Tc of the 
thicker NbN islands in the bilayer. 


3.4. Bi 2 Se 3 — NbN bilayers 


We now present the main part of this study, namely, 
the properties of a highly resistive 10 nm 812863 on 3 
nm NbN bilayer, and this after an extensive preparatory 
part on the properties of films and bilayers pertinent 
to the present investigation. It turned out that the 
results of the R and MR in Figs. 5 and 6 of a similar 
bilayer (on half a wafer) are not very interesting, since 
the resulting MR above Tc is very similar to that of a 
bare 10 nm 812863 film. This is a direct result from 
the fact that the links between superconducting islands 
of the 3 nm reference NbN film were still too strong, 
as is obvious from the almost full transition to zero 
resistance at low temperatures. We therefore had to 
weaken the inter-grain links even more. Since the use 
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FIG. 11: (Color online) Resistance versus temperature of 
a pristine 10 nm Bi 2 Se 3 on 3 nm NbN on FS bilayer under 
different magnetic fields and gate voltages. The inset shows 
conductance spectra of this highly resistive bilayer at 1.86 K, 
0 T and -100 Vg. 

of 2.5 nm thick NbN films on FS yielded very resistive 
films, with R over a Mfl at low temperatures, with a 
lot of noise and a very weak signature of a transition 
at 2-3 K, we decided to stick with the 3 nm thick films 
in the bilayers. We thus prepared and characterized 7 
wafers with the 10 nm Bi 2 Se 3 on 3 nm NbN bilayers, 
and found out that their maximum resistance (at low 
temperatures) varied between about 1 and 10 Ml. This 
depended mostly on the base pressure used before the 
deposition process (1.5-3x10“^ Torr), but probably also 
on the surface quality of the FS wafers. In any case, the 
most interesting results were obtained with the wafers 
of highest resistance, and here we present the results on 
the most resistive one, which also has the weakest links 
as we originally planned to have. 

Fig. 11 shows the resistance versus temperature of 
this wafer with the 10 nm Bi 2 Se 3 on 3 nm NbN on 
FS bilayer. The temperature dependence of R clearly 
shows insulating behavior down to 2.5 K where a small 
decrease of resistance is observed. Similar behavior 
was observed also in 5 nm thin cuprate films |25j . We 
show data at four locations on the wafer which exhibits 
quite a spread, and also the data of contact C6 at 0 



FIG. 12: (Golor online) Magnetoresistance versus temper¬ 
ature of the C4 contact of the pristine bilayer of Fig. 11, 
obtained for 0 and 1 T fields and under gate voltages of 0 and 
±100 Vg. All other contacts show similar behavior. The inset 
shows the full temperature range while the main panel is a 
zoom in on the low temperature regime. The inset also shows 
the MR of the reference 10 nm thick Bi 2 Se 3 film of Fig. 6. 

and 1 T magnetic field under 0 Vg. Clearly, the small 
decrease of R below 2.5 K disappears under 1 T, and 
we thus conclude that it originates in superconductivity 
of the NbN islands which are very weakly linked. The 
inset shows conductance spectra of all these contacts 
at 1.86 K which clearly exhibit tunneling behavior 
between the superconducting islands. The gap-like 
feature appearing at about 70 meV, is a result of several 
weak-links connected in series between the voltage 
contacts. Assuming an energy gap A of 2 meV for the 
NbN islands [2S] yields 35 such weak-links in series. The 
dips in the conductance spectra of the inset are due to 
heating effects as a result of reaching the critical current 
of the NbN islands m- 

The main results of this study are presented in Fig. 
12 and its inset, where the MR of the bilayer of Fig. 
11 for 0 and 1 T and under different gate voltages 
shows a complex, non monotonous peaks structure. For 
comparison, the corresponding monotonous MR of the 
reference 10 nm Bi 2 Se 3 film of Fig. 6 is also shown. 
Below 2.5 K, the MR is similar to the results in the 
inset to Fig. 5, but with a 3-5 times smaller signal at 
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the maximum MR at ~1.9 K. The later reflects the fact 
that the weak-links here are much weaker than in Fig. 
5. There is also no visible effect of the gate voltage 
on the MR results in this regime. Above 2.5 K, the 
MR data becomes very sensitive to the gate voltage, 
and up to 7-8 K exhibit kind of anti-phase behavior 
versus temperature for ±100 Vg where the peak and 
dip at 3.5 K reverse roles at 6-7 K. The MR under 0 
Vg is somewhere in between these two at 3.5 K, but 
then becomes closer to the -100 Vg data at 6-7 K. At 
higher temperatures the MR data becomes even more 
complex, but clearly it goes down at 25-30 K, where it 
tends to follow the reference Bi^Se^ film. It should be 
stressed here that the MR behavior of Fig. 12 for the C4 
contact, was observed also for the other contacts on this 
wafer (C6 and C8), while C7 had an additional strong 
oscillatory behavior above 10 K as shown in Fig. 13 for 
0 Vg. This oscillatory behavior was washed out under 
±100 Vg. The inset to Fig. 13 is a zoom-in on weaker 
oscillations or plateaus in the MR of C7 in the range 
of 3-8 K. Such plateaus and knees, although weaker, 
appear also in the other contacts. As far as we know, 
there is no theory that predicts oscillations of MR versus 
temperature in a 2D topological superconductor. 

To find out whether these oscillations and plateaus 
originate in magnetic field effects, we measured R versus 
H for all the contacts on the wafer of Fig. 11 at 1.85 
K. The results are depicted in Fig. 14, and they are 
quite similar to those of the inset to Fig. 6, besides the 
difference in overall resistance and aging effects. We 
have also measured R versus H under 0 Vg on this wafer 
(when it was 1.5 month old) at higher temperatures. 
A background of a polynomial of order 2 fit to the 
data was subtracted from the data, and the resulting 
resistance variation Ai? versus H was obtained. The 
normalized AR(H)/R(0) results at ^9.5 K are plotted 
in Fig. 15 versus magnetic field. Presuming that the 
observed oscillations originate in flux quantization in 


T(K) 



FIG. 13: (Color online) MR versus temperature at OVg 
of the C7 contact on the wafer of Fig. 11. The oscilla¬ 
tions were smeared out under ±100 Vg. The inset is a zoom 
in on the data at low temperatures. Only this C7 contact 
showed the pronounced oscillations versus temperature, al¬ 
though plateaus and less pronounced knees were observed also 
in the other contacts as seen in the inset and marked by the 
arrows. 



H(T) 


FIG. 14: (Color online) Resistance versus magnetic field of 
the bilayer of Fig. 11 at 1.85 K and 0 Vg. Although this bi¬ 
layer was kept in dry air for one month prior to this measure¬ 
ment, it deteriorated (aged) as can be seen from the increase 
of its resistance values at 0 T by about 4 kfl as compared to 
Fig. 11. 


nano-holes in the bilayer, we obtain from the period 
A7J ~ 0.6 T a hole diameter of 66 nm, which agrees 
quite well with the AFM image and line profile of 
Fig. 1. This result suggests that superconductivity 
in the bilayer persists at least up to ~9.5 K. Similar 
oscillations were observed also in the other contacts at 
3.7 and 5.3 K, but they disappeared at 15.1 K in the 
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0.6-1 10nm Bi2Se3/3nm NbN 



Tir^ = (b /AH ~ 3.4x10'''® ==> r{hole) - 33 nm 

-0.6 J^^^-n---T-—- ^^^^-. 

0 1 2 3 4 5 6 

H(T) 

FIG. 15: (Color online) Oscillations of normalized R minns 
background [AR(H)/R(0)] of the C7 contact of Fig. 13 versns 
field at 9.48 K and 0 Vg. 

noise level of the measurements. Thus Tc of the NbN 
islands in this bilayer should be between 9.5 and 15.1 K. 
This is in agreement with the the MR onset in Fig. 10 
under 4 T (12.2 K), or the temperatures of maximum 
resistance of Fig. 9 (12-13 K), although in that case 
the NbN film was 4 nm thick and not 3 nm thick as in 
the present bilayer. We conclude from these findings 
that Tc of the superconducting islands in our bilayer can 
reach 12-13 K, and that if fluctuations are taken into 
account, superconductivity could persist in our bilayers 
up to ~15 K. Thus, it should be possible to observe 
flux-flow effects of vortices up to this temperature, but 
this however, does not explain the non monotonous MR 
features observed in Figs. 12 and 13. Above 15 K, any 
MR effect different from that due to the reference 10 nm 
Bi 2 Se^ film, should have a different origin. 

Next, we discuss possible explanations for the complex 
behavior of the MR results of Fig. 12. In principle, 
the present bilayer as drawn schematically in Fig. 
4, comprises of six distinct superconducting regions, 
marked in this figure by I to VI. The NbN islands (I) 
and the inter-islands regions (II) are the obvious ones, 
but proximity induced superconductivity occurs also 
in the Bi 2 Se^ layer above these regions, in particular. 


at the bottom and top surfaces of this doped topolog¬ 
ical layer where the helical edge currents flow. This 
yields four more superconductive regions in the cap 
R* 2 «S'e 3 layer, two above the islands (III and V) and 
two above the inter-islands regions (IV and VI). Due 
to proximity effect and distances from the supercon¬ 
ducting NbN islands (I), we can safely assume that 
rc(I)> rc(III)> Te(II)> Tc(IV). We shall assume first 
that the MR results observed in Fig. 12 originates in 
vortex physics, and discuss them in terms of pinning in 
the different superconducting regimes I-VI. Then, for the 
remaining unidentified features, we shall discuss possible 
contributions to the MR from zero energy Majorana 
modes and unconventional superconductivity. 

In the context of vortex physics, the MR peaks seen 
in Fig. 12 as a function of temperature originate in flux- 
flow, and are due to two competing phenomena. One is 
the increased vortex generation on cooling down due to 
newly formed proximity induced superconducting areas 
in the bilayer which eventually saturates, and the other 
is increased pinning with decreasing temperature. The 
result of these effects is an MR peak versus temperature 
below which all vortex motion stops and the MR goes 
to zero. Such MR peaks have been clearly observed 
in similar cuprate bilayers [28], as well as in the inset 
to Fig. 5 here. Since the bottom Bi 2 Se 3 surface of 
regions III and IV are very close to regions I and II, the 
proximity induced superconductivity in them is stronger 
than in the more distant regions V and VI. We shall 
first discuss the MR versus T result of Fig. 12 under 
0 Vg. The MR of this data is most always lower than 
that of the reference Bi 2 Sez film. This looks like MR 
suppression in the bilayer, but since there are zero MR 
regimes in the bilayer data (at 6-7 and 18-20 K), it looks 
more like an indication of strongly interacting layers 
in the bilayer. This interaction seems to play a lesser 
role only above 25 K where the two curves coincide, 
apparently because no superconductivity affects the MR 
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results at this high temperature range. We attribute the 
double MR peak between 7-15 K to flux flow, vortex 
generation and pinning in regions I and III. We note 
here that the above analysis was facilitated by the fact 
that this double peak is fully separated from the other 
features of the MR curve. Using similar arguments, 
we can attribute the broad MR knee between 2.5 and 
6.5 K to pinning in regions II, and the large MR peak 
below 2.5 K to pinning effects in region IV. The fact 
that the MR peak below 2.5 K is so strong might be 
due to enhance proximity induced superconductivity by 
the surface currents of the topological layer in region IV, 
or to additional contributions from other effects, to be 
discussed below. 

We shall now analyze the MR results of Fig. 12 
under different gate voltages and try to explain the 
two prominent peaks at 3.5 and 6 K under -100 Vg 
and -1-100 Vg, respectively. We note that transport, 
and therefore also the MR data, is strongly dependent 
on the conductance of regions II and IV, since they 
comprise the weak-links for current flow in the bilayer. 
As a result, they will be more sensitive to gating than 
region I. Under -100 Vg, the depletion layer in Fig. 4 
in between the NbN islands (region II), is positively 
charged or simply electron depleted. This makes region 
II, (as well as the adjacent region IV just above it), 
less conducting (or more insulating), thus lowering its 
Tc onset to about 4 K compared to its value without 
gating (onset at ~6.5 K). The MR peak at 3.5 K under 
-100 Vg is consistent with this scenario, and therefore 
seems to originate in region II. The large MR peak at 
1.9 K seems to be unaffected by either of the present 
gate voltages of ±100 Vg, and this again might originate 
in topological effects in region IV. Under ±100 Vg, the 
depletion layer is electron rich (more conducting), thus 
its proximity induced T^. in region II is higher than the 
4 K obtain before under -100 Vg. The MR peak at 6 K 
under ±100 Vg can thus originate also in region II. Fig. 


12 shows extra MR features under gating besides the 
two prominent peaks at 3.5 and 6 K. These include the 
broad hump under ±100 Vg between 10 and 25 K, and 
the broad peak between 15 and 30 K under -100 Vg. As 
discussed before, superconductivity dies off at about 15 
K, and any proximity induced will obviously be lower. 
Therefore, if one wishes to invoke vortex motion as the 
origin of these extra MR features, one has to assume the 
existence of a pseudogap phase for the present bilayers. 
Signature of such a phase was found in point contact 
conductance spectra measurements of copper doped 
Cuo. 2 Bi 2 Se 3 single crystals with Tc ~3 K, where a 
depletion of density of states at low bias was observed up 
to ^20 K |2S]- Considering that in the present bilayers 
the NbN islands have a much higher ~12-13 K, the 
conjectured pseudogap phase could easily reach 25-30 
K. Hence, to summarize, the extra MR features in Fig. 
12 above 15 K could originate in flux flow, provided a 
pseudogap phase exists in our system. In the cuprates, 
where the pseudogap phase is well established, such 
vortices were detected in thermoelectric measurements 
under a magnetic field (the Nernst effect) (SHI- 

Finally, we focus on the intensity or magnitude of the 
MR peaks in Fig. 12, and on what might contribute 
extra strength to the dominant peak at 1.9 K. We note 
that the volume fraction of the superconducting NbN 
islands with the highest T^. values is apparently very 
small, as can be inferred from the very small resistive 
transition at 2.5 K in Fig. 11. Previously, in the 
vortex pinning context, we attributed the MR double 
peak under 0 Vg at 7-15 K to regions I and III, the 
knee between 2.5 and 6.5 K to region II, and the large 
MR peak below 2.5 K to region IV. It was hard to 
reconcile how region IV, with its presumably weakest 
superconductivity, leads to such a large MR peak. Thus, 
if the enhanced MR peak at 1.9 K is still due to pinning 
effects, we are led to the conclusion that region IV must 
have enhanced superconductivity just as well. Such en- 


12 


hancement might originate in a longer normal coherence 
length in the Bi 2 Se 3 at low temperatures, which 
is possibly enhanced further by hybridization of the 
helical surface currents of regions IV and VI. We know 
of no theory that predicts such effects, and therefore 
this interpretation is only a hypothesis at the present 
time. An alternative scenario for the interpretation of 
our results is that the MR knee between 2.5 and 6.5 K 
is due to pinning effects in both region II and region 
IV. This would leave the dominant MR peak at 1.9 
K unaccounted for, and other effects as for its origin 
should be explored. But this time, vortices in the bilayer 
are fully pinned and presumably form a vortex lattice. 
Vortices in general, in a topological superconductor, are 
predicted to host zero energy Majorana modes in the 
vortex cores [2]. Thus, in a dense periodic vortex lattice, 
interactions between these modes lead to the creation 
of new energy bands in the band gap of this system, 
some of which could be quite flat [31]. Whether or 
not these Majorana bands contribute to the MR of the 
bilayer is unknown at the present time. Rigorous MR 
calculation versus temperature are needed in order to 
test if this scenario actually occurs in the present bilayers. 

4. CONCLUSIONS 

A comprehensive study of ultra thin bilayers com¬ 
prising of weakly connected s-wave superconducting 
islands and a continuous topological insulator cap-layer 
reveal interesting gate and temperature dependent 
magnetoresistance features. These features could be 
explained as originating in vortex physics, where prox¬ 
imity induced superconductivity and pinning effects in 
the different regions of the bilayer play a major role. 
In particular, the present results are consistent with 
enhanced superconductivity in the surface region of the 
topological layer in contact with the superconductor, 
where the helical edge currents flow. Signature of a 
pseudogap was found, and a possible contribution to the 


data from Majorana bands was discussed. 
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